function determ(array,norder)
      implicit double precision (a-h,o-z)
      dimension array(10,10)
10    determ=1.
11    do 50 k=1,norder
      if (array(k,k)) 41,21,41
21    do 23 j=k,norder
      if (array(k,j)) 31,23,31
23    continue
      determ=0.
      goto 60
31    do 34 i=k,norder
      save=array(i,j)
      array(i,j)=array(i,k)
34    array(i,k)=save
      determ=-determ
41    determ=determ*array(k,k)
      if (k-norder) 43,50,50
43    k1=k+1
      do 46 i=k1,norder
      do 46 j=k1,norder
46    array(i,j)=array(i,j)-array(i,k)*array(k,j)/array(k,k)
50    continue
60    return
      end